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Abstract 

A thin shell of light-like dust with its own gravitational field is studied in 
the special case of spherical symmetry. The action functional for this system 
due to Louko, Whiting, and Friedman is reduced to Kuchaf form: the new 
variables are embeddings, their conjugate momenta, and Dirac observables. 
The concepts of background manifold and covariant gauge fixing, that underlie 
these variables, are reformulated in a way that implies the uniqueness and 
gauge invariance of the background manifold. The reduced dynamics describes 
motion on this background manifold. 
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1 Introduction 



The phenomenon of gravitational collapse leads to serious problems in the classi- 
cal theory of gravity. The structure of the resulting singularities contradict the 
foundations of the theory such as the equivalence principle. Thus, the existence of 
singularity theorems |J may constitute a strong motivation to address the quanti- 
sation of the gravitational field, in the hope that such a framework can avoid the 
occurrence of singularities. 

A prominent feature of the classical collapse is the existence of horizons which 
appear sooner than the singularity. Such horizons not only imply that the singularity 
is inevitable (which is, roughly, the content of the singularity theorems), they also 
seem to prevent any object or information from leaving the region of collapse and 
from coming back to the asymptotic region. It is the existence of horizons that 
makes gravity so different and the problem of collapse so difficult. On the other 
hand, the problem of gravitational collapse is a very special one. For its solution, 
a complete quantum theory of gravitation may be as little needed as the complete 
quantum electrodynamics was needed for the first calculations of atomic spectra. 

Motivated by these ideas, we consider the quantum theory of a spherically sym- 
metric thin shell and its gravitational field. This is, in fact, a quite popular system. 
For example, it was used to study the motion of domain walls in the early Universe 
0, of black- hole evaporation H, of quantum black holes H, and many others. In 
H, || and @, gravitational collapse of such a thin shell in its own gravitational field 
has been studied. The result was analogous to what is known about the s-mode of 
the Coulomb problem. Two aspects of this result were surprising. First, for low- 
mass shells, there were stationary states with Sommerfeld spectrum and scattering 
states with the wave packets describing the shell bouncing and re-expanding. The 
evolution was unitary. Second, there was an analogue to the critical charge in the 
relativistic quantum mechanics of atoms. The role of charge was played by the rest 
mass of the shell (not to be confused with its total energy) and the critical value of 
the rest mass was about one Planck mass. As in the case of relativistic atoms, the 
quantum-mechanical description breaks down for supercritical "charges" . 

To understand these results was difficult. Even the simple scattering of the sub- 
critical shells admitted several different interpretations. There were two problems. 
First, the radial coordinate of the shell, which served as the argument of the wave 
function, did not possess the status of a quantum observable. The Coulomb-like 
potential prevented one from constructing a position operator similar to the Newton- 
Wigner operator. Second, the model was completely reduced to the physical degrees 
of freedom, which in this case is just the radius of the shell. However, the value of 
the radius is not as informative in a black-hole spacetime as it is, for example, in 
Minkowski spacetime: points with the same value of radial coordinate can lie in 
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different asymptotically flat regions. These regions are separated by horizons. The 
tools that were at one's disposal in H, H] and |7j did not allow to decide whether 
the shell created a horizon and then, consequently, re-expanded behind this horizon 
into a different asymptotically flat section, or whether it did not create any horizon 
and re-expanded into the same region from which it collapsed. 

In ||, two remedies have been proposed. The first is to work with a null (light- 
like) shell. The classical dynamics of such a shell is equivalent to that of free photons 
on flat two-dimensional spacetime (the "charge" is zero). For such a system, there 
is a well-defined position operator ||. Moreover, it admits a simple description of 
its asymptotic states, unlike Coulomb scattering. 

The second idea is that the equation of motion for r(t), which has been obtained 
from Einstein's equations, must in fact result from a reduction to true degrees of 
freedom of an action that contains the shell as well as the gravitational field. Indeed, 
in the spherically-symmetric case, the gravitational degrees of freedom consist only 
of the gauge and the dependent ones. It seems that the reduction has been performed 
in such a way that the information about the geometry of spacetime has been lost. 
We shall, therefore, perform the reduction explicitly in a careful way. This is the 
main purpose of this paper. 

In general, the reduction procedure consists of two steps: the choice of gauge 
and the solution of constraints. There exists a particular form of the gravitational 
action that is effectively reduced, but which still contains some information about 
the geometry of spacetime: the so-called Kuchaf decomposition [K|, [ft!]] . Kuchaf 
variables are neatly separated into pure gauge ones (so-called embeddings), depen- 
dent ones that are conjugate to the embeddings, and physical degrees of freedom. 
Some progress in understanding Kuchaf decomposition has been achieved recently 



12] , general existence of the decomposition has been shown, and the crucial 



role of gauge choice in it has been recognised. The nature of gauge choice in quan- 
tum gravity has also been elucidated. Two important notions have been introduced: 
background manifold and covariant gauge fixing. As a matter of fact, the present 
paper is the first practical application of these concepts. It deals with the classical 
canonical analysis of the null-dust shell. Its main result is the explicit construction 
of the Kuchaf decomposition. This then serves as the starting point for the quantum 
theory of the shell, which will be presented in a separate paper. 

The plan of the paper is as follows. Sect. 2 explains the notions of background 
manifold and covariant gauge fixing in a new and clear form. This enables us to 
show the gauge invariance, the uniqueness and some additional structures of the 
background manifold, which will be necessary for the interpretation of the shell 
quantum mechanics. Some important points of ]TB[ are then summarised in the 
new language. Sect. 3 describes the solutions of Einstein's equations containing the 
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shell. After fixing the gauge, these solutions are reformulated as a set of parameter- 
dependent metric fields and shell trajectories on the background manifold. The 
parameters distinguish the physically different solutions and will play the role of the 
physical variables in the Kuchar decomposition. The representative metric fields 
and shell trajectories will be used to define the transformation from the ADM to 
Kuchar variables on the constraint hyper- surface. This transformation is performed 
in Sect. 4. Starting point is a (non-reduced) Hamiltonian action principle [14| for 
the spherically-symmetric shell and its gravitational field. In Sect. 4, an extension 
of the results to a whole neighbourhood of the constraint hyper-surface is performed 



using the methods and theorems of [13]. The final action, the variables in it and 
some discussion can be found in Sect. 5. 



2 Background manifold 

and covariant gauge fixing 

In this section we introduce the two basic notions of background manifold and 
covariant gauge fixing, restricting ourselves, for the sake of simplicity, to vacuum 
general relativity. The language is slightly different from that used in so that 
we can prove more results; we also summarise some points that are important for 
the paper. 

Let Ai be a four-manifold that admits a Lorentzian metric field g such that (Ai, g) 
is a globally hyperbolic spacetime. Dynamically maximal solutions of Einstein's 
equations are always of this form [|1J. Then according to a theorem of Geroch 
fl6fl , Ai = £ x R, where £ is an initial data manifold. Topological sectors of 
general relativity are uniquely associated with different three- manifolds S. Ai is 
called background manifold. In this way, each topological sector determines a unique 
background manifold At. 

All diffeomorphisms of Ai form the group DifLM; this is considered as the "gauge 
group" of general relativity. Observe that the group depends on the topological 
sector chosen, i.e., on Ai, and that the manifold structure of Ai itself is gauge 
invariant. Single points of Ai are, however, not gauge invariant, since they are 
being pushed around by DifLM. In some important cases, not DifLM but some of 
its subgroups play the role of gauge group. For example, in the case of asymptotically 
flat space-times, only those diffeomorphisms are considered as gauge transformations 
that become sufficiently quickly trivial at infinity. In general, in such cases, Ai is 
equipped with some gauge-invariant structure in addition to the naked manifold 
one. 

Let ip G DifLM and let g be a Lorentzian metric on Ai. Then the inverse pull- 
back associated with ip maps g into another metric g', g' = ((f~ 1 )*g- In this way, 
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the group DifLM acts on the space RienxM of all (suitably restricted) Lorentzian 
metrics on Ai. The action is not transitive and so there are non-trivial orbits of 
the group in RiemA'l. Such orbits are called geometries and the quotient space 
RiemA'l /DifLM is the space of geometries. Let tt : RiemA'l i— > RiemA'l /DiffAl be 
the natural projection for the quotient. One can equip the space of geometries with 
some additional structure, e.g., a topology, starting from a structure of RiemA'l and 
using the projection. 

Suppose that we manage, at least for some open set U C RiemA'l /DiffAl, to 
specify a section a. This is a map, 

a : RiemA'l /DiffAl h-> RiemAl 

such that 7roo" = id. The meaning of such a section is that a particular representative 
metric on AI is chosen for each geometry in U. This is exactly what has been called 
covariant gauge fixing in |13[| . Clearly, the transformation between two covariant 
gauge fixings is not a single diffeomorphism, but an element of the Bergmann-Komar 
group 
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Given a covariant gauge fixing a on U, one can use it to construct a map from 
RiemAl /Di&M x Emb(S, Ai), where Emb(E, Ad) is the space of embeddings of the 
initial data surface £ into Ai, to the ADM phase space of general relativity. The 
construction has been described in and it goes, roughly, as follows. Let 7 G U 
and let (7(7) be the representative metric on Ai. Let X : £ 1— > Ai be an embed- 
ding. Then (7(7) determines the first, qui{x), and the second, Kki(x), fundamental 
forms of the surface A(£) in the spacetime (Ai,a(j)). The corresponding point 
(qki{x), ir kl (x)^j of the ADM phase space can be obtained in the well known way 
from q kl (x) and K kl (x). 

In |[L3|| , this transformation has been restricted to give only the points of the 
constraint surface V of the ADM phase space; moreover, only those points of V have 
been selected, where the evolved space-times do not admit any isometry. Then the 
map from RiemAl /Diff Ai x Emb(£, Ai) to V has been shown to be invertible and 
extensible to a neighbourhood of RiemAf/DifLM x Emb(S, Ai) in the larger space 
RiemAl /DifLM x T*Emb(S, Ai), which has then been mapped to a neighbourhood 
of T in the ADM phase space. Next, the Darboux-Weinstein theorem has been 
employed to prove some nice symplectic properties of the map. These properties 
then make the map to a general transformation of the ADM to the Kuchaf variables. 

This procedure will here be applied to the model of spherically-symmetric thin 
gravitating shell in the subsequent sections. We shall find in the next section the 
set of representative solutions for Einstein's equations for each physically distinct 



situation of the shell because, as has been shown in [13], this part of the section 



a suffices completely to construct the above map to the constraint surface of our 
model. 
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3 Einstein dynamics of the shell 



Any spherically-symmetric solution of Einstein's equations with a thin null shell as 
the source has a simple structure. Inside the shell, the spacetime is flat; outside the 
shell, it is isometric to a part of the Schwarzschild spacetime of mass M. The two 
geometries must be stuck together along a spherically-symmetric null hyper-surface 
so that the points with the same values of the radial coordinate R coincide. 

All physically distinct solutions can be labeled by three parameters: rj e {— 1, +1}, 
distinguishing between the outgoing {rj = +1) and in-going {rj = —1) null surfaces; 
the asymptotic time of the surface, i.e., the retarded time u = T — R e (— oo, oo) 
for 7] = +1, and the advanced time v = T + R e {—oo, oo) for rj = — 1; and the mass 
M e (0, oo). An in-going shell creates a black- hole (event) horizon at R = 2M and 
ends up in the singularity at R = 0. The outgoing shell starts from the singularity 
at R = and emerges from a white-hole (particle) horizon at it! = 2M. 

We can write down the metric in the case rj — 1 with the help of retarded 
Eddington-Finkelstein coordinates U, R, d and (p. U = u is the trajectory of the 
shell, U > u is a part of Minkowski spacetime, 

ds 2 = -dU 2 - 2dUdR + R 2 dn 2 , (1) 

and U < u is a part of Schwarzschild spacetime, 

ds 2 = - (l - ™) dU 2 - 2dUdR + R 2 dfl 2 . (2) 

Similarly, for rj = — 1, the advanced Eddington-Finkelstein coordinates are V , R, d 
and ip, and V = v is the shell. Inside the shell, V < v, 

ds 2 = -dV 2 + 2dVdR + R 2 dn 2 , (3) 

and outside the shell, V > v, 

ds 2 = - (l - dV 2 + 2dVdR + R 2 dfl 2 . (4) 

Let us denote the spacetime given by the triple of parameters 77, M and w by 
{rj, M, w), where w = u for rj = 1 and w — v for 77 = — 1. 

We observe that the two space-times {r],M,wi) and {rj,M, w 2 ) are isometric, 
the isometry sending the point {U ,R,d,(p) into ([/ + u>2 — u>i, -R, tp). Hence, the 
geometries of the solutions that differ only in the value of the parameter w are equal. 
Yet, the physical situations they represent are different; this is similar to the motion 
of a free mass point in Minkowski spacetime. For each two different trajectories, 
there is a Poincare transformation that sends the first into the second. Still, the two 
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motions are physically different because they look differently from one fixed inertial 
frame. For the shell, instead of an inertial frame, we imagine that there is a fixed 
asymptotic family of observers. The group of these isometries is a symmetry group 
rather than a gauge group. It can (and will) be employed to define a time evolution. 

Another interesting isometry is the map T : (77, M, w\) 1— > (—77, M, w 2 ) defined for 
77 = +1 and arbitrary w-y and w 2 by the Eddigton-Finkelstein coordinates as follows: 

TiU^R^A^i) = (V 2 ,R 2 ,# 2 ,<p 2 ) , 

where 

V 2 — —Ul + Wl + W 2 , R 2 —Rl, &2 = '&l, ( f>2 = ( fl- 

For rj = — 1, we just take the inverse of the above so that T 2 = id. T can be viewed 
as a time reversal symmetry. 

The Eddington-Finkelstein coordinates may be nicely adapted to the symmetry 
and may simplify the metric, but they do not define a covariant gauge fixing. Indeed, 
the identification of the points (Ui, R\, $1, </?i) of the solution (+1, Mi, w±) with the 
points (V2, R 2) $2, f2) of (— l,M 2 ,w 2 ) satisfying the relations V 2 = U\, R 2 = R±, 
$2 = $1 and ip 2 = ipi will invert the time orientation of the asymptotic observers, 
which is to stay gauge invariant. We need, however, a covariant gauge fixing if we 
are to transform the action to the Kuchar form. The rest of this section will be 
devoted to a choice of gauge that will be convenient for this problem. 

To start with, we have to specify the background manifold. Our model com- 
prises only the spherically-symmetric part of general relativity with the shell. We 
shall, therefore, admit only spherically-symmetric initial surfaces X and only that 
subgroup of DifiM (where M := ExR), the elements of which commute with 
the rotations and are trivial at infinity. Let p, $ and <p be coordinates on E that 
are adapted to the symmetry and p e [0, 00), where p = is the regular centre of 
symmetry; we assume that E is smooth at this centre. The shell is at p = r, and 
the infinity at p = 00. 

The coordinates d and tp are ignorable coordinates; in the action, we can in- 
tegrate over them so that they disappear and the effective initial manifold E is 
one-dimensional, diffeomorphic to R + , and the effective background manifold Ai is 
two-dimensional, R + x R. Our restricted gauge group induces an effective gauge 
group, Diffo )00 .M, on A4; it only contains diffeomorphisms that preserve the central 
boundary as well as, pointwise, the infinity. 

Let us choose coordinates U and V on M. that satisfy the following boundary 
conditions at the gauge- invariant boundaries of M.: At the regular centre inside the 
shell, 

U = V , 
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at Z~, U = — oo and V G (—00, 00), at X + , V = 00 and [/ G (—00, 00), and at i°, 
U = —00 and V = 00. Otherwise, {7 and V are arbitrary. 

Using these coordinates U and V, one can define the representative metric (see 
Sect. 2) by conditions on its components with respect to U and V . We shall choose 
them as follows. 

I. U and V are double- null coordinates so that the representative line element 
takes the form 



ds 2 = -A(U, V)dUdV + R 2 (U, V){d$ 2 + sin 2 tfV) • 



(5) 



2. The representative metric is continuous at the shell. 

3. For the outgoing shells, U is the retarded time determined by the representative 
metric at V = 00. Analogously, for the in-going shells, V is the advanced time 
at U — —00. 

Such a metric is uniquely defined for any physical situation given by the values 
of the parameters 77, M, and w. This can be shown as follows. 

Consider first the case rj = +1. The Eddington-Finkelstein coordinate U satis- 
fies already the conditions for U, so we need only to find the function V. In the 
Minkowski part, U > u, of the solution, the boundary conditions at the centre lead 
uniquely to: 



A = 1 



R 



V-U 



(6) 



In the Schwarzschild part, U < u, of the solution, V is an advanced null coordi- 
nate, so it must be some function, V = X(M, u, V), for each fixed M and u, of the 
advanced Eddington-Finkelstein coordinate V, which is defined by 



V := U + 2R + 4Mln 



R 



2M 



- 1 



The function X is uniquely determined by the boundary condition at the shell, 
requiring that V be continuous: 

X(M,u,V\ u=u ) = (U + 2R) u=u , 



or, 



X (^M,u,u + 2R + 4M\n 
To solve this equation, we define 



R 



2M 



= u + 2R. 



x := U + 2R + 4Mln 



R 

2M 
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calculate R in terms of M, u, and x, and substitute the result into the right-hand 
side: 

X(M,u,x) = u + 2R(M,u,x). (7) 
A straightforward calculation yields 



R(M,u,x) = 2M«(exp( 



x — u 



\ r v 4M 

where k is the well-known Kruskal function defined by its inverse, 

K- 1 (y) = (y-l) e y , (9) 
and R > 2M was used. Eqs. (|7|) and (g) yield 

((A_ l)exp( ™)). (10) 

A similar calculation for R < 2M leads to the same result. From this, it is easy to 
calculate R, if we observe that 

R +lnf— -lUlnfV 1 ' R 



2M \2M J \ \2M 
Then Eq. (|) implies that 

— ((^-0-(^r))- 

This relation defines the desired transformation from (U, R, ip) to (U, V, tp). 

As the last step, we calculate the metric for U < u in the new coordinates. First, 
we differentiate the function R. The derivative of ([]) determines the derivative of 



which holds for any /. Then, 



dR ' 



2 K (/ + )e K (/+) 
with 



y-u \ (V-U\ JTT V-u fV-U. n . 

1 ex P 777 + - T77 - ex P - 77T - ^ 



4M J L \ AM J AM 1 \ AM 



Now, (|IT]) implies that 

«(/+) R 
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and (§) that / = \K(f) — lj expf/c(/)). This leads to the relation 



/+ x _ 2M 



K{f+)e<f+) R 



and thus 



1 / 2M\ Jrr 1 V-u V -U\ 
dR= — 1 — dll H -— — ttt — — exp — — — dV . 

2 V R J 2K(f + )e<M 4M F \ AM 1 

Substituting this into the metric @ results, finally, in 

R = 2Mn{f + ), A= fr ] . f V ~ r U exp ( —^-) , (15) 
w+y ' n{f + )e<f+) AM \ AM 1 



where /+ is defined by (0), cf. (fTT|) . With these expressions, it is easy to verify that 
A and R are continuous at the shell, as required. We note that these expressions 
contain u as well as M, which become conjugate variables in the canonical formalism. 
This makes the transition to the embedding variables non-trivial, and one must first 
look for this transformation on the constraint surface. 

In the case of in-going shells (rj — —1) a completely analogous procedure yields, 
for V < v, again (|6]), and for V > v, 

i2 = 2Af «(/_), A= - I r . , V ~^ exp (YjlE-\ (16) 
w ; K(/_)e K (/~) AM F V AM J ' v ; 

where 

, fv-U \ (V-U 
f- '■= —rvr- - 1 exp 



AM J V AM 
These expressions result from (|I5D by the substitution V — u — > v — U. 

As the result of the gauge fixing, the set of solutions (77, M, w) can be written as 
a set of (r), M, w)-dependent metric fields (|) and a set of shell trajectories on a fixed 
background manifold M.. Here, the corresponding functions A and R have the form 

A( V , M, w; U, V) , R{r], M, w; U, V) , (17) 

and the trajectory of the shell on the background manifold is simply U = u for 
rj = +1 and V — v for rj = — 1. 

A key property of the background manifold is that it possesses a unique asymp- 
totic region with I~ defined by U — > —00 and I + by V —>■ +00. As the shell cannot 
escape the background manifold, its reappearance at an asymptotic region must be 
interpreted as the reappearance at the asymptotic region of M.. In this way, the 
background manifold is a tool to solve the problem of where the shell reappears. 
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4 Transformation to embedding variables 



4.1 Canonical formalism 

The form of the canonical theory that is based on the embedding rather than ADM- 



type variables has been studied and advocated by Kuchaf. In the recent paper [13 
a large step forward in this field has been achieved. The embedding variables have 
been associated with background manifolds and gauge fixings similar to what has 
been done in the previous section. The existence of this transformation has been 
shown in the general case. 

The resulting formalism inspires hopes that some unpleasant features of the ADM 
variables can be removed. First, the ADM variables lead to singular points in the 
physical configuration space (super-space |L8|, as well as at the constraint surface 



corresponding to spaces or space-times with symmetries. Second, the symmetry 
of the ADM theory itself is, on one hand, too large, containing all infinitesimal 
surface deformations, including also those transformations that do not result from 
diffeomorphisms. On the other hand, it is too small because only infinitesimal 
surface deformations and not finite group elements can act on the whole phase 
space. The constraint surface that has been constructed in has, however, the 
form of a fibre bundle, which is a manifold (all points are regular), and the fibre 
group of this bundle is the diffeomorphism group of the background manifold, so it 
acts on the whole bundle. 

As a Hamiltonian action principle that implies the dynamics of our system, we 
take the action Eq. (2.6) of |14[] (see also ||). Let us briefly summarise the relevant 



formulae. The spherically symmetric metric is written in the form: 
ds 2 = -N 2 dr 2 + A 2 {dp + N p dr) 2 + R 2 dn 2 , 

and the shell is described by its radial coordinate p = r. The action reads 
S = Jdr pr + J dp(P A A + P R R - H 

and the Hamiltonian is 

H = NH + N p H p + N^E^ , 

where := lim^^ N p (p), is the ADM mass (see fTJ), N and N p are the 
lapse and shift functions, 7i and 7i p are the constraints, 

AP 2 P\Pr RR" RR'A' R' 2 A m> , 
H p = P R R' - PIA - p5(p - r) , (20) 
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and the prime (dot) denotes the derivative with respect to p (r). 

The main topic of this paper is to transform the variables in the action Sq. This 
transformation will be split into two steps. The first step is a transformation of the 
canonical coordinates r, p, A, Pa, R, and Pr at the constraint surface V that is 
defined by the constraints (|19D and (p0|). The new coordinates are u and p u = —M 
for r] — +1, v and p v = —M for 77 = —1, and the so-called embedding variables U(p) 
and V(p). 

The second step is an extension of the functions u, v, p u , p v , U(p), Pu(p), V{p), 
and Pv(p) ou t of the constraint surface, where the functions u, v, p u , p v , U(p), 
and V(p) are defined by the above transformation, and Pu(p), Pv{p) by Pu(p)\r — 
Pv{p)\r = 0. The extension must satisfy the condition that the functions form a 
canonical chart in a neighbourhood of T. A proof that such extension exists in 
general has been given in [T3I . 



4.2 Transformation functions at the constraint surface 

The constraint surface contains only points of the phase space that correspond to 
initial data for solutions of Einstein's equations. Hence, we can assume that the 
metric (0) is a spherically-symmetric solution with a shell, and so the functions 
A(r], M, w] U, V) and R(r), M, w; U, V) are those written down in the previous section, 
Eqs. (H), (|15|) and (0). According to Sect. 2, if such a metric is given, then, for each 
embedding, a unique first and second fundamental form can be calculated from it, 
and so the map from the embeddings to the ADM variables qki{x) and n kl (x) can 
be constructed. 

A very important point is to specify the family of embeddings that will be used 
throughout the paper. The embeddings are given by 

U=U(p), V = V(p). 

These functions have to satisfy several conditions. 

1. As S is spacelike, U and V are null and increasing towards the future, we must 
have U' < and V > everywhere. 

2. At the regular centre, the four-metric is flat and the three-metric is to be 
smooth. This implies U'(0) = — V(0) in addition to the condition U(0) = 
V(0). This follows from T'(U(0),V(0)) = and means that S must run 
parallel to T = const, in order to avoid conical singularities. 

3. At infinity, the four-metric is the Schwarzschild metric. We require that the 
embedding approaches the Schwarzschild-time-constant surfaces T = const, 
and that p becomes the Schwarzschild curvature coordinate R asymptotically. 
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More precisely, the behaviour of the Schwarzschild coordinates T and R along 
each embedding U(p), V(p) must satisfy 

T(p) = T^ + Oip- 1 ), (21) 
R(p) = p + 0(p- 1 ). (22) 

The asymptotic coordinate is a gauge-invariant quantity and it possesses 
the status of an observable. 

4. At the shell (p = r) we require the functions U(p) and V(p) to be C°°. In fact, 
as the four-metric is continuous in the coordinates U and V, but not smooth, 
only the C 1 -part of this condition is gauge invariant. Jumps in all higher 
derivatives are gauge dependent, but the condition will simplify equations. 

We suppose further that there is a whole foliation of the solution space-times. 
Any foliation can be considered as a one-parameter family of embeddings: 

U = U(r,p), V = V(t, P ), 

the parameter being r. The metric (EH reads, in terms of the coordinates r, p, # 
and <p: 

ds 2 = -AUVdr 2 - A{UV + VU'jdrdp - AU'V'dp 2 + R 2 dQ 2 . 

From this metric, we can read off the values of the variables A, R, N p and A^ 
immediately: 

A = y/-A(o,U,V)U'V', R = R(o,U,V) } (23) 
where o symbolises the observables (w and M, respectively), and 

2U'V 2U'V 

The expression UV — VU' is the Jacobian of the transformation from r and p to U 
and V, and we assume it to be positive. 

To calculate the gravitational momenta, we can use the canonical equations that 
follow when the action So is varied with respect to Pa and Pr: 

P A = -|(P-AW), 

Pa = -|(P-iV p P')-f (A-(iV p A)') • 
Substituting for R, A, A^ and N p gives 

A - (N p A)' = ^t UV '~Z U ' (-AuU' 2 V + A V U'V' 2 - AU"V + AU'V" 
UV - VU' 

R — N P R' = 2WVI (RuU* - R V V) , 
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so that, finally, 



11 m (RuW - R V V) , (24) 



y/-AU'V 



RA Utt> RA Vt „ RU" RV" 



P R = R V V - R V Y + -rftf - + 1W -^- (25) 

Here, the indices U and V denote the partial derivatives with respect to U and V. 

Eqs. (0), (^4[), and (p5[ ) are the transformation equations expressing the variables 
A, i?, Pa and Pr in terms of the new variables at the constraint surface. The 
functions A and R are given by fllBD , and ([U]). 

We now turn to the remaining variables 77, r and p. We let rj unchanged; in fact, 
we shall consider the action as two different actions, one for each value of rj. The 
variable r is related to our new variables u, v, M, U(p) and V(p) in a different way 
for each value of rj. If rj = +1, then r is determined by the equation U(r) = u. This 
is an equation with exactly one solution if u satisfies the condition tt < {7(0) because 
U(p) is a monotonous function with the range (—00, {7(0)). For the differentials of 
the variables U(p), r and u, we obtain the relation: 

du — dU(r) , , 

*=-i^r- (26) 

Similarly, if 77 = —1, then r is defined by V(r) = v for v > V(0), and the relation 
between the differentials takes the form: 

du - dV(v) 

V'(r) ' 1 j 



The variable p does not seem to be determined completely in |U] because equa- 



tion (2.5a) of Ref. Hli| , which is the only equation that could serve this purpose, 



does not make sense in the limit m — > of null shells. However, the constraint 
equations lead to some expressions for p; these determine p only at the constraint 
surface, but this is, in fact, all we need. Let us, therefore, turn to the constraint 
equations. 



4.3 The constraints 

The constraint functions flT9|) and (|20|) contain finite parts, which are obtained for 
p 7^ r and in the limits p — > r±, and ^-function parts. The 5-function parts can be 
rewritten as equations for finite quantities, if one collects all terms with 5-function 
and sets the coefficient equal to zero. 

From the boundary conditions at the shell and Eqs. (pBp, ( |2~4] ) and (P5p, it follows 
that the functions A(p) and R(p) are continuous, whereas A'(p), R'(p), -Pa(p) and 
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Pr{p) jump across the shell, as the metric is not smooth. This implies in turn that 
the ^-function part of the constraints is equivalent to 

p = - V R{R'} , (28) 
P = -A[P A ] , (29) 

where the symbol [g] := g+ — g~ denotes the jump of the quantity g across the shell. 
Let us calculate the jumps. We have 

[R'} = [Ru}U' + [R V }V . 

For 77 = +1, we have to use ([]) and ([15]) and to replace the limits p — > r± by 
U — > m±. We obtain immediately from that 

Ru- = -\, Rv- = \- 



Differentiating (O) with the help of formulae ( |T2"D and flT3] ) leads to, for U < u, 

V-u pvri V-u 

p _ p _ 4M 4M 

-TV/+ — — : — r-r-, — -77— s , -Kv 



f pyti y_u 

■/ + p = 4M eX P 



Eq. (D and = R/2M imply for the limits that 

lim K(f+) = , 

so we have 

1 2M 1 

2 V — u 2 

Hence, 

[^1 = ^-, [i?y]=0. 

Similarly, for 77 = —1, 

[ifc]=0, [Ry] = --^L. 

There is also the relation 

K\u=u — — - — , n\v=v — — 2 — 

and so (ESI) yields: 



77 = +1 : p = -MU'(r) , (30) 
r] = -1 : p = -MV\r) . (31) 
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For Pa? (0) implies 

A [P A ] = R[Ru}U' - R[R V ]V 



and so (p9|) gives the same result as (pq). 

Let us return to the finite part of the constraints (|T9| ) and (p0|). If we substitute 
the above expressions for A, P, Pa, and Pr, we obtain, after some lengthy but 
straightforward calculation, for each p/r: 

H = ^{ARR UV + ARuR v + A)U'V+ 

R ' ' " ' — ^-/2 i ^ f AD A D M/'2 



— (AP^ - A^P^f/' 2 + ^(Ai2 vv - A v P y )V 



and 



ft p = -^(AP^ - AuRu)U' 2 + -{AR VV - A V R V )V' 2 



^AR uu -A u R u )U' 2 + ^[ 

If 7i and 7i p are zero for any embedding outside the shell, that is for all possible U' 
and V, the coefficients of U'V, U' 2 and V' 2 must themselves vanish: 

ARRuv + ARtiRv + A = 0, (32) 
ARuu-AuRu = 0, (33) 
AR VV -A V R V = 0. (34) 

These three equations are equivalent to the full set of Einstein equations for any 
metric of the form (|5]). Thus, our functions A and P have to satisfy these equations. 
This is immediately clear for (|B|) which gives A and P inside the shell. A more 
tedious calculation verifies the validity of (^D, (|33|), and ( |34"|) also outside the shell, 
where A and P are given by fllED and ( |T6f) . 



4.4 Transformation of the Liouville form 



As it has been explained at the end of Sec. |4.1| , the transformation of the action 
(|18|) to the new variables will be performed in two steps. The first step is restricted 
to the constraint surface and forms the content of the present section. 
At the constraint surface, H = A^P^ and the action ([TJ) becomes 

S \r = Jdr pr - N^E^ + J dp{P h k + P R R) 



According to the discussion given in [2~0|, the ADM boundary term N^E^ in the 
action can, after parametrisation at the infinity, be written as E^T^ and can be 
considered as a part of a modified Liouville form. Let us denote this form by G: 



6 = J dp (P A dA + P R dR) + pdr - E^dT a 



(35) 
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As a result, the transformation of the action is nothing but the transformation of 
the Liouville form 0. 

We expect that the terms remaining after the transformation do not depend on 
any embeddings, because the pull-back of the symplectic form to the constraint 
surface is degenerated exactly in the direction of the gauge variables U(p) and V(p). 
As we shall see, the constraint surface F consists of two components, T + and r~, 
T + containing all outgoing and T~ all in-going shells. We split this form into three 
terms for the in-going and outgoing part, respectively, 

6| r + = + | r + + 6~| r + + pdr, 

where 

/oo 
dp (P A dA + P B dR) - MdT^ 

because — M at the constraint surface, and 



- | r+ = j T dp (P\dA + P R dR) 
Jo 



and similar expressions for 0|r-- Let us first transform the part 0|r+ of the Liouville 
form and make the ansatz 

/•OO 

0+|r+ = J dp§ - MdToo , (36) 

with 

d = (fdU + gdV + hide?)' + dip , (37) 

where we have denoted the observables u and M collectively by o l (i = 1,2). This 
has to be compared with the corresponding part of (j35|), where the substitutions are 
made from (p3|), 

dA Au JTT a. Ay AM a. A A < _l dU ' a. dV ' (V* 

T = ^ du+ ^ dv+ 2A d0+ W> + W<' (38) 

dR = R v dU + R v dV + R i do i , (39) 



and Pa and Pr are given by (p4|) and (|25|). 

It turns out to be convenient to make for the functions in (13?T) the ansatz 



/ = ^(-y)+mv,j), (40) 

9 = ^\a{~^+G{UK<t), (41) 

hi = + H^V^) t (42) 

(p = RRuU' -RR V V -^(RuU' + R v V')\n^-^J (43) 

-FU' - GV + (j)(U, V, o*) . (44) 
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The functions F, G,Hi,(f) are then determined through comparison with the coeffi- 
cients of dll, dV, dU', dV, and do 1 . This leads to the equations 

Fy-Gu = ^(2AR UV - A V R V - A V R V ) , (45) 

Hiu — Fi — -^{2AR iV - AiRu - AyR,) , (46) 

Htv-Gi = —(2ARi V — AiR v — A v Ri) , (47) 
= 0. (48) 

We next calculate the right-hand side of these equations by using the explicit ex- 
pressions for A and R found in Sect. 2. Outside the shell, these are the expressions 
(15). It is convenient to introduce the abbreviations 

. V-u U-u 



AM 1 AM 

One then has 



(49) 



h e b-a 

A = , R = 2Mk , (50) 

where k is a function of /+ = (6 — l)e b ~ a , see (13). The following identity turns out 
to be useful: 

„b—a ,„ 1 

(51) 



e K 6-1 

After some lengthy, but straightforward calculations one finds 

Fv -°« = -8^1) ' (52) 
H uU = F u-^^Yy (53) 

H -v = G u + -£f - . , (54) 



86(6 


-1) ' 


K - 


- 1 


86(6 


-1) ' 


1 


K — 


2 ~~ 


2(6- 


K 


- 1 



Hmu — Fm — — — 7777 — , (55) 

k 2 - 6 2 

Hmv = ° M + 26(6^1) a + 26(6 = 1) • (56) 

The freedom in the choice of solution to Eq. (fE>D— (fI7D enable us to set F = 0. From 
the first equation one then gets 

L dU W^) = - 46(6 - 1) ' (5?) 

where we have chosen the boundary condition that G = for U = u, i.e., at the 
shell, and calculated the integral by the substitution x = k. One recognises from 
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(pwj) that, at the shell, Gm = and G u = —1/86. With this result for G, one can 
integrate Eqs. (|53|)-(l5"6|) for Hi and choose the integration constants such that 

H u — —G , (58) 
H M = ~(U-u)-4bG, (59) 

having Hi = 0, i — 1, 2, at the shell. This then yields for the Liouville form outside 
the shell pulled back to the constraint surface 

0+| r+ = (fdU + gdV + hidd)^ + d dp <pj + dvy\ p=r - MdT^ . (60) 



The fourth term on the right-hand side of (|60|) is a total derivative and will be 
omitted, since it does not contribute to the dynamics. Eqs. (|26"|), (|30|) and ([H]) lead 
to 

pdr = -M(du - dU) . 
Analogously, one finds for the part 0~|r+ inside the shell 



e- 



r+ 



(kdU + ldV)\ r + d (J* dp tf?j - dr^j\ p=r , (61) 



with 



i(j = RRuU' - RR V V - ^(RuU' + R v V')hi^~A . (64) 

Compared to /, g, hi, <p, there are no terms analogous to G, F, and Hi, since the 
classical solutions (|6|) inside the shell lead to a vanishing right-hand side of (f45|)- 
((47]). Because of the boundary condition U'(0) = — V'(0), the functions k and / 
vanish at the centre. The third term on the right-hand side of ( |6l~D is again a total 
derivative and will be neglected. 

One has therefore only potential contributions at the shell and at infinity. We 
shall consider first the contribution from the shell. Since there F = G = H u = 
H M = 0, one has to calculate 

(kdU + ldV)\ p=r - drip\ p=r - (fdU + gdV + h u du - h M dM)\ p=r 
+dvif\ p=r - M(du - dU) . (65) 

Using (§) and (|i~5D , one arrives at the following jump conditions at the shell: 

[RRu] = M, [RR V ] = 0, [RR U ] = —M, [RR M ] = . (66) 
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Taking these into account, one recognises that all terms on the dust shell cancel. 
As we shall now demonstrate, the only non-vanishing terms are originating from 
infinity. 



ei 



r+ 



In 



lim 

p— >oo 

' -u r 



HM + FdU + GdV + 



^dU + ^dV + ™±dd 



MdT n 



(67) 



where the function F — 0, and G, H u and H M are given by fl57|), (5S|), and 
respectively The limit ( |67D is determined by the boundary conditions 3 of Sec. 4.1, 
cf. (f?Tp and (0). 

Eqs. (^) and (^) determine the expansions of U(p) and V^(p) uniquely. Indeed, 
for 7] — +1, U near the space-like infinity coincides with the Edington-Finkelstein 
retarded coordinate (see Sect. 3) and so is given in terms of T and R by 

R 



U = T - R-2M\n 



2M 



Then, 



U{p) 



-p-2M\n ^)+T 0O + 0(p- 1 ). 



(6£ 



The presence of the logarithmic term is due to the long range of the gravitational 
potential. Thus, the first diverging term is universal, the second depends on the 
observable M, and the asymptotic coordinate of the embedding appears only at 
the third position. 

The asymptotic expansion of the function V(p) can be determined from (|Il]). We 
have first to get rid of k: 



(V -u- AM) exp 



V 

AM 



2(R- 2M) exp 



( U 



R 



\AM 2M 



(69) 



Then we substitute the expansions (22) and (68) into the right-hand side of (^ 

V 



(V — u — AM) exp ■ 



AM 

(p - 2M + CK/T 1 )) exp 



AM 



p-2M\n^ + T oo )+0(p- 1 ) 



(70) 



Let us remove the singular part in the exponent by setting 

V(p) = p - 2M In + + V x {p) 



8M 



Eq. ( |70"D then becomes 

(1 -2Mp- 1 )expO(p- 1 ) 



1 - 2Mp~ l In + (T^ - u - AM)p~ 1 + V^p' 1 



exp 



Vi(p) 
AM 



(71) 
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Taking the limits p — > oo of both sides of (|7T|), we obtain 

p^oo I p ) 1 AM 



It follows that 

lim Vi(p) = 



Eq. (|7lD implies then also that 

Hence, the expansion of the function V"(p) has the form: 

V{p) = p-2M\n^ + T oo + ; (72) 

the asymptotic coordinate appears again only at the third position, and this is 
the reason why the expansion must be carried so far. 

Now, the expansion of all functions contained in is a straightforward matter. 
For G, we obtain from (|57): 

M AR 2 - (V - u) 2 



G 



A (V -u) 2 -AM(V -u)' 
Then, Eqs. (J2|), (H), and (|7|) give 

G = --M - AM 2 p~ 1 In ^ + M(2T ao - 2u - ZM)p~ l + ofjT 1 ) , 

where o(p _1 ) is defined by the property lim^oo po(p~ l ) = 0, and 

H u = -G = -M + AM 2 p~ 1 In - M(2T 00 -2u- 3M)p~ 1 + 0(p- 2 ) 

as well as 

H M = ^ p+ 7 -Mln^ + ^(6 + 41n2) - - u) + O^ 1 ) . 

Eq. (|72|) can be used to calculate dV: 

^ = d(-2Min^) + ar. + o(^) . 

Then we obtain: 

HM + FdU + GW = M(dT 00 - dw) + rfZ + O ^— , 
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where 

Z = -- Af -u) + ^Mp + ^M 2 In p + i(4 - 13 In 2)M 2 - ^M 2 In M. 
Similarly, 

fln(^)=2M + 0(p- 2 ). 

The derivatives ifo, Ry, R u and -Rm can be expanded if we calculate them from 
flljj) using the identity (pl|): 

\R-1M 1R-2M V-u 

Ru — — — — - — , -R 



2 i? 2 i? V-M-4M 

R-2M 1 
i4 = -2M- 



and 



This gives: 



R V-u-iM 



R R - 2M ( 1 V - [/ - u 



~ /l/T 2 D \ AMI 4M 



M i? V4Ml-^2- y-«-4M 



R In ( J (ifodl/ + iZyciV + ifcdo') = O ' lu / ! 



2 v y/ 7 VP/ 

Collecting all terms, we finally have: 

6| r+ = -Mdu + dZ + O (^-J . (73) 

The exact form dZ can be omitted because it has no influence on the symplectic 
form and the equations of motion. 

The final result of this subsection can be formulated as follows. The constraint 
surface T consists of two components: r + for the outgoing shells (rj = +1), and 
for the in-going shells (77 = —1). On r + , we have the coordinates M, u, U(p) and 
V(p), and the pull-back of the Liouville form to T + is 

9| r+ = -Mdu. 

Thus, it is independent of U(p) and V(p), as expected. 

In a completely analogous manner, the following result can be derived for the 
7] = — 1 case: 

9| r _ = -Mdv , 

and our coordinates on T~ are M, v, U(p) and V(p). 

These results also show that the two Dirac observables — M and u (or — M and 
v) form a conjugate pair. Indeed, the Poisson algebra of Dirac observables is well- 
defined by the (degenerate) pull-back of the symplectic form to the constraint sur- 
face. 
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4.5 Extension to a neighbourhood of the constraint surface 



In the previous subsection, the constraint surface pull-back of the Liouville form 
has been transformed to the Kuchaf coordinates: the embeddings U(p) and V(p) 
that represent pure gauges, and p u and u (or p v and v) that are Dirac observables. 
The next task is to extend these coordinates to a neighbourhood of the constraint 
surface r = T + U r~. 



In [13], the proof has been given that such an extension exists if there are no 



points with additional symmetry at T. For the action (18), the space-time solutions 
are all spherically symmetric, but none of them exhibits any other symmetry, be it 
discrete or continuous. We can reformulate the result of |T3|] in a way suitable for 
our purposes as follows. There is a neighbourhood T f± of in the phase space, 
and functions 

U,Pu,V,P v ,u,p u , (74) 



in T /+ and 



in T' such that, at T z 



U,Pu,V,P v ,v,p v , 



(75) 



0. 



Pu(p) = P v (p) 

and U(p), V(p), p u and u (or p v and v) coincide with our coordinates there. The 
functions ([F4]) and ( [75|) form canonical charts in r /=t . The transformation between 
the old variables 

R,P R ,A,P A ,r,p, (76) 

and the new ones Eq. ([74]) or ([75| ) is smooth and invertible in r . At the constraint 
surface r*, the transformation is given by Eqs. (p3|)-(l25|), (IB]), (|T6|), U(r) = u, 
V(r) = v, fl3~0|) and (|31|). Outside the constraint surface, only the existence of the 
transformation has been shown so we do not know its form. 

Using this result, we can write the transformed action in r as follows, 



dr 



Pn ii-\- I dp(P u (p)U(p) + 



Pv{p)V{p) - N u {p)P u {p) - N v {p)P v {p 



and 



S' 



p v v+ dp(p u {p)U{p) + 



dr 



P v {p)V{p) - N u (p)P u (p) - N v {p)P v {p) 



The two actions can be considered as the reduced form of just one action. Con- 
sider the case 77 = +1 first. The dynamical trajectory of the shell is given by 
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the relation u(r) = const, whereas v(t) is arbitrary, depending on the choice of 
the parameter r (the only restriction is that v{r) is an increasing function). This 
information can be obtained from the extended action: 



s; xt = I dr 



p u ti + p v v - n u p v + j dp (Pu(p)U (p)- 
P v (p)V(p) - N u (p)P u (p) - N v {p)P v {pi 



where n u is a Lagrange multiplier , v (r) a pure gauge and ^-dependent. Similarly 
for f] = —1: 

/oo poo . 

dr p u u + p v v - n v p u + / dp [P u (p)U(p) + 
-oo l_ JO 

P v (p)V(p) - N u (p)P u (p) - N v {p)P v {p)) . 
One can set in S^ xt , as p u = —M < 0: 

n u = np u 

and, similarly, in S~ xt , 

n v = np v . 

Then, clearly, S£ xt and S~ xt are obtained by reducing the following action 

/oo roo . 

dr p u ii + p v v - np u p v + dp [P u (p)U(p)+ 
-oo L Jo 

Py{p)V{p) - N u (p)P u (p) - N v {p)P v {p)) . 

Indeed, the case r] = +1 (77 = — 1) is obtained from the solution p v = {p u = 0) of 
the constraint 

PuPv = . 

The relation between the total energy M and the two momenta p u and p v can 
now be written as follows: 

M = -p u -p v . 

5 Conclusions 

We have demonstrated that there is a transformation of variables bringing the action 
(P^) to the simple form of the so-called Kuchaf decomposition 

p p poc 

S = Jdr {p u ii + p v v - np u p v ) + Jdrj Q dp{P v U + P V V - H) , (77) 
where H = N u Pjj + N v Py] n, N u (p) and N v (p) are the new Lagrange multipliers. 
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The dependence of the new variables Pjj and Py on the old ones is not known. 
This dependence would be needed for calculation of the spacetime geometry associ- 
ated with any solution given in terms of the new variables. We know, however, that 
the new constraint equations, Pu{p) — Pv{p) — 0, are mathematically equivalent to 
the old constraints, Eqs. ( [19]) and (|2~0l) . One can, therefore, use the old constraints to 
calculate the geometry from the true degrees of freedom along the hypersurfaces of 
some foliation. The fact that two spacetimes obtained by this method using differ- 
ent foliations are isometric is guaranteed by the closure of the algebra of constraints 
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The new phase space has non-trivial boundaries: 

p u < 0, p v < , (78) 

^>0, (79) 

p v = , U E (— oo, u) , V > u — Ap u K ( — exp — ) , (80) 



and 

p u = , V E (v, oo) , U < v + Ap v K — exp 



4p u 
V-v" 

4p v , 



The boundaries defined by inequalities ([79| ) -(|81| ) are due to the singularity. 

The two dynamical systems defined by the actions (|1^) and ( |77| ) are equivalent: 
each maximal dynamical trajectory of the first, if transformed to the new variables, 
give a maximal dynamical trajectory of the second and vice versa. 

The variables u, v, p u and p v span the effective phase space of the shell. They 
contain all true degrees of freedom of the system. One can observe that the cor- 
responding part of the action (|77|) coincides with the action for free motion of a 
zero-rest-mass spherically-symmetric (light) shell in flat spacetime if one replaces 
the inequality (|79|) by 

—u + v 

> . 

2 

Such a dynamics is complete if the singularity at the value zero of the radius of 
the shell, (— u + v)/2, can be considered as a harmless caustic so that the light can 
re-expand after passing through it. It might, therefore, seem also possible to extend 
the phase space of the gravitating shell in the same way so that the in-going and 
the out-going sectors are merged together into one bouncing solution. 

However, such a formal extension of the dynamics ( |77] ) is not adequate. The 
physical meaning of any solution written in terms of new variables (|74|) or (75) 



is given by measurable quantities of geometrical or physical nature such as the 
curvature of spacetime or the density of matter. These observables must be expressed 
as functions on the phase space. They can of course be transformed between the 
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phase spaces of the two systems ( |T8| ) and (|77|). They cannot be left out from any 
complete description of a system, though they are often included only tacitly: an 
action alone does not define a system. This holds just as well for the action ( |T8"D as 
for (0). 

Let us consider these observables. The expression for the stress-energy tensor of 
the shell written down in [JTJJ implies that the density of matter diverges at r = 0; 
this corresponds to (— u +v)/2 — in terms of the new variables. Eqs. and 
(|j~6l) can be used to show that the curvature of the solution spacetime diverges at 
the boundary defined by Eq. (|80|) for p v = and by Eq. fl8T|) for p u = 0. It follows 
that the observable quantities at and near the "caustic" are badly singular and that 
there is no sensible extension of the dynamics defined by action ( |T?D to it, let alone 
through it. This confirms the more or less obvious fact that no measurable property 
(such as the singularity) can be changed by a transformation of variables. 

The action for the null dust shell is now written in a form which can be taken as 
the starting point for quantisation. Surprisingly, it will turn out that the quantum 



theory is, in fact, singularity-free. This will be done in a separate paper 22 
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